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1.  Summary 


This  report  discusses  the  research  conducted  by  the  U.S.  Anny  Research  Laboratory,  Vehicle 
Technology  Directorate  (ARL-VTD)  on  advanced  suspension  control.  ARL-VTD  research  on 
advanced  suspension  systems  should  result  in  a  reduction  of  ground  vehicle  chassis  vibration 
while  maintaining  tire  contact  with  the  road  surface.  The  purpose  of  this  research  is  to  reduce 
vibration-induced  fatigue  to  the  Warfighter,  as  well  as  to  improve  the  gun  shooting  precision. 


2.  Introduction 


Passive  suspension  systems  for  ground  vehicles  are  optimized  for  passenger  comfort  or 
road/terrain  handling.  Quite  often  these  two  goals  contradict  each  other.  Modern  suspension 
systems  may  include  control  systems  that  use  sensors  and  actuators  to  improve  both  perfonnance 
metrics  at  the  expense  of  increased  complexity  and  cost. 

The  two  main  groups  of  advanced  suspension  systems  are  divided  into  fully-active  and  semi¬ 
active  actuators.  Semi-active  suspension  systems  vary  spring  and  damper  properties  to  store  or 
remove  energy  from  the  system  according  to  the  applied  control  algorithm.  Most  commonly, 
this  is  achieved  through  variable  damping  components:  variable  orifice  valves  or  magneto- 
rheological  (MR)  fluids.  Fully-active  suspension  systems  use  actuators  to  add  and  remove 
energy  from  the  system.  This  is  achieved  using  hydraulic,  pneumatic,  or  electromechanical 
actuators.  This  report  focuses  on  fully-active  actuator  control  algorithms. 

The  baseline  control  algorithm  for  this  research  is  the  Linear  Quadratic  Regulator  (LQR).  This 
algorithm  has  been  extended  to  incorporate  future  road  profile  disturbance  data  through  Preview 
Control.  Another  algorithm,  Generalized  Predictive  Control,  was  evaluated  with  formulations 
for  Implicit  Disturbances  and  Explicit  Disturbances.  Generalized  Predictive  Control  with 
Preview  Control  was  not  considered  in  this  completed  research  work.  This  will  be  considered  in 
future  research. 


3.  Modeling  and  Simulation 


Two  primary  vehicle  models  were  developed  in  Mathworks  Matlab-Simulink  to  gain  insight  into 
the  behavior  of  the  vehicle’s  suspension  subject  to  various  road  conditions.  These  models  were 
idealized  representations  of  a  quarter  of  the  vehicle.  The  models  were  implemented  in 
Mathworks  Matlab/Simulink  for  simulation. 


1 


The  Quarter  Car,  as  shown  in  figure  1,  is  modeled  as  25%  of  the  chassis  mass  riding  on  a 
suspension  system  consisting  of  a  spring,  damper,  and  actuator.  This  is  connected  to  the  wheel, 
which  is  modeled  as  a  lumped  mass  with  a  very  stiff  spring  and  very  soft  damper.  Road 
disturbance  is  introduced  through  the  wheel  suspension  system.  The  control  actuator  is  treated  as 
a  perfect  force  generator;  thus,  there  are  no  dynamics  associated  with  the  actuator  model. 


Figure  1.  Quarter-car  suspension  model. 


In  the  representative  model  shown  above,  the  lumped  masses  each  have  a  single  degree-of- 
freedom  -  vertical  motion  of  the  mass  centroid;  therefore,  this  is  a  two-degree-of-freedom  (2 
DOF)  system.  This  system  is  then  modeled  as  a  system  of  second  order  differential  equations. 
In  this  linearized  model,  two  assumptions  must  be  noted: 

(1)  The  tire  always  remains  in  contact  with  the  ground,  and 


(2)  There  are  no  stroke  limits  for  the  suspension  components. 
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[i]  =  A  [x]  +  Bc[u]  +  Bd [d] 


where 


[x]  is  the  state  vector  comprised  of:  zv,  zs,  zu,  and  zs, 

zv  is  the  height  of  the  tire  centroid  above  a  starting  reference  zero, 
zs  is  the  height  of  the  chassis  centroid  above  a  starting  reference  zero, 
z'u  is  the  vertical  velocity  of  the  tire, 
zs  is  the  vertical  velocity  of  the  vehicle  chassis 


2 


[u]  is  the  control  vector  comprised  of  each  control  actuator 

[d]  is  the  disturbance  vector  comprised  of:  d,  and  d, 

Subscript  c  indicates  the  control  component  of  the  input  matrix,  and 

Subscript  d  indicates  the  disturbance  component  of  the  input  matrix 

The  particular  parameters  chosen  for  each  component  in  the  Quarter  Car  modeled  are  listed  in 
tables  1  and  2.  These  are  taken  from  two  papers  discussing  recent  active  control  research.  The 
two  sets  of  vehicle  parameters  will  be  referred  to  as  Light  Vehicle  (2)  and  Heavy  Vehicle  (5). 
These  models  are  used  to  reproduce  data  for  the  control  systems  in  this  report.  This  provides  a 
common  reference  point  for  comparison  of  the  control  algorithms;  these  do  not  necessarily 
represent  parameters  for  a  specific  vehicle. 

Table  1.  Light  Vehicle  Parameters  used  in  research  by  El  Madany,  et  al.,  (2). 


Vehicle  Parameters  -  Light  Vehicle 

ms 

1155.6  kg 

total  chassis  (sprung)  mass 

mu 

28.58  kg 

tire  (unsprung)  mass 

Ks 

19960  N/m 

suspension  spring  rate 

Cs 

816.46  Ns/m 

suspension  damping 

Ku 

155900  N/m 

tire  spring  rate 

Cu 

“0^ 

tire  damping 

1st  mode 

1.237  Hz 

chassis  natural  frequency 

2nd  mode 

12.15  Hz 

tire  natural  frequency 

stiffness  and  damping  chosen  to  provide  damping  ratio  of  0. 17 

Table  2.  Heavy  Vehicle  parameters  used  in  algorithm  research  by  van  der  Aa,  et  al.,  (3). 


Vehicle  Parameters  -  Heavy  Vehicle 

ms 

34600  kg 

total  chassis  (sprung)  mass 

mu 

1350  kg 

tire  (unsprung)  mass 

Ks 

440000  N/m 

suspension  spring  rate 

Cs 

43100  Ns/m 

suspension  damping 

Ku 

6500000  N/m 

tire  spring  rate 

Cu 

“0^ 

tire  damping 

1st  mode 

1.055  Hz 

chassis  natural  frequency 

2nd  mode 

10.95  Hz 

tire  natural  frequency 

4.  Controller  Methodology 


Several  different  controllers  were  considered  and  employed  on  the  suspension  models  developed. 
The  performance  of  the  controllers  was  evaluated  against  each  other  with  the  passive  suspension 
model  as  the  baseline.  These  algorithms  are  multi-input/multi-output  controllers.  Thus,  the 
control  command,  [ut],  can  be  a  vector  where  each  element  corresponds  to  a  value  of  a  specific 
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actuator  at  time  t.  The  same  is  applicable  for  the  measured  response,  \yt\,  and  external 
disturbances,  [dt\. 

All  of  the  algorithms  discussed  in  this  research  are  classified  as  optimal  control  algorithms. 
Optimal  control  algorithms  calculate  a  control  command  that  best  satisfy  a  specified  objective 
function.  The  objective  function  is  often  a  linear  combination  of  weighted  state-vectors, 
measured  responses,  and  control  commands  integrated.  The  minimum  of  the  objective  function, 
found  by  setting  the  derivative  of  the  objective  function  to  zero,  corresponds  to  the  optimum 
command.  The  weighting  matrices  serve  two  purposes.  First,  they  are  used  to  normalize  and  non- 
dimensionalize  the  input/output  data.  Thus,  large  accelerations  and  small  displacements  can  be 
used  to  provide  a  dimensionless  objective  function.  The  second  purpose  of  the  weighting 
matrices  is  that  they  allow  relative  performance  penalties.  For  example,  high  accelerations  in  the 
crew  cabin  should  be  penalized  more  than  the  same  level  of  acceleration  in  the  cargo  bay. 

These  algorithms  are  formulated  as  regulators.  Regulator  Control  derives  the  control  command 
required  to  drive  the  feedback  inputs  to  zero.  Conversely,  Tracking  Control  drives  the  feedback 
inputs  to  follow  a  specified  path.  Deadbeat  Regulator  Control  derives  the  control  command 
such  that  for  a  discrete-time  model  the  closed  loop  poles  are  zero  in  the  complex  plane;  thus,  the 
closed  loop  system  is  critically  damped. 

Two  other  concepts  are  applicable  to  algorithms  in  this  research.  Preview  Control  or  Look-Ahead 
Control  measures  the  external  disturbance  before  it  affects  the  vehicle  and  incorporates  this 
information  into  the  control  algorithms.  LQR  with  Preview  Control  is  discussed  below.  System 
Identification  is  a  process  where  an  algorithm  identifies  the  vehicle  model  based  on  a  history  of 
measured  inputs  and  outputs.  Various  control  laws  can  then  be  calculated  and  applied  to  the 
vehicle.  Generalize  Predictive  Control  algorithms  perform  a  system  identification  on  the  fly  to 
calculate  the  control  laws. 


4.1  LQR 

Linear  Quadratic  Regulator  is  the  most  basic  of  the  optimal  control  algorithms.  This  control 
algorithm  is  a  full-state  feedback  controller.  Given  a  perfect  model  of  the  vehicle  dynamics,  the 
optimal  control  command  can  be  calculated  and  applied.  Full-state  feedback  requires  that  every 
state  in  the  model  can  be  measured  and  fed  into  the  control  law. 


Recall  the  state  space  model  of  the  vehicle  dynamics  from  the  previous  section: 

[x]  =  A  [x]  +  Bc[u]  +  Bd[d] 

For  LQR,  the  control  command  is  calculated  as  a  matrix  gain  multiplied  by  the  state-vector. 
Thus,  the  control  command,  u,  is  calculated  by 

[u]  =  -K[x], 


where 
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K  is  the  LQR  feedback  matrix  gain,  and 
/.r/is  the  state  vector. 

The  feedback  matrix  gain  is  obtained  by  building  a  quadratic  objective  function,  J,  and  solving 
for  the  minimum: 

n  co 

/  =  I  (xTQx  +  uTRu  +  2  xT  Nu)dt 
Jo 

where 

Q  is  a  symmetric  positive  semi-definite  matrix  that  assigns  weights  to  the  state  vector, 

R  is  a  symmetric  positive  definite  matrix  that  assigns  weights  on  the  control  commands, 

N  is  a  symmetric  positive  semi- definite  matrix  that  assigns  weights  for  the  cross-coupled 
terms.  N  is  often  assumed  to  be  null. 

The  superscript  T  indicates  the  vector  or  matrix  transpose 
The  minimum  is  obtained  by  taking  the  derivative  of  the  function  and  setting  it  equal  to  zero. 

0  =  -J-  =  xtQx  +  utRu  +  2xtNu 
at 

The  equation  devolves  into  an  Algebraic  Riccati  Equation  (ARE),  and  the  LQR  feedback  gain 
matrix  is  calculated  from 

K  =  R-'BlP, 

where  P  is  the  solution  to  the  ARE. 

As  implemented  in  this  research,  LQR  feedback  control  is  formulated  for  discrete  time  domain. 

4.2  LQR  with  Preview  Control 

Linear  Quadratic  Regulator  with  Preview  Control  extends  the  LQR  with  an  additional  control 
command  component  based  on  future  disturbance  data.  Preview  Control  was  first  published  by 
Bender  (1 ),  and  is  further  developed  in  papers  by  El  Madany,  et  al.,  (2)  and  van  der  Aa,  et  ah, 

(2). 

The  closed  loop  state  matrix  is  used  to  calculate  the  optimal  command  for  each  disturbance  in  the 
future.  The  resulting  commands  are  integrated  to  provide  a  weighted  augmentation  to  the  LQR 
closed  loop  control  command,  as  shown  below: 

[U]  =  -K[X]  -  R~'BTC  Review  ^A-BCK)^ 

where  the  preview  horizon,  tPreview,  is  the  time  that  corresponds  to  the  furthest  valid  future 
disturbance  data  point. 
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The  number  of  samples  of  future  disturbance  can  be  varied  based  on  the  capability  of  the  sensor 
and  vehicle  ground  speed.  Again,  for  this  research,  this  control  algorithm  was  formulated  in 
discrete  time  domain.  Performance  of  the  control  system  for  the  chosen  number  of  Preview 
samples  was  compared. 


4.3  GPC  with  Implicit  Disturbance 

Generalized  Predictive  Controller  (GPC),  as  detailed  in  a  paper  by  Kvatemik,  et  al.,  (4)  and  the 
text  by  Juang  (5),  is  a  linear,  time-invariant,  multi-input/multi-output  controller  that  generates  an 
Autoregressive  with  Exogenenous  (ARX)  model.  This  means  that  GPC  builds  a  time-domain 
based  model  from  time -history  data  of  measured  responses  and  control  commands.  This  ARX 
model  is  used  to  derive  a  control  law  that  implements  an  optimal  deadbeat  regulator.  GPC  is 
capable  of  performing  this  system  identification  and  control  law  generating  the  model  in  near- 
real-time.  The  system  identification  portion  requires  that  the  controller  inject  frequency-rich 
control  commands  into  the  system,  typically  band-limited  white  noise  (random  signal  with  a  flat 
power  spectral  density).  The  issued  commands  and  resulting  system  responses  are  used  to  build 
the  ARX  model.  External  disturbances  are  not  directly  measured  and  are  implicitly  accounted  for 
in  the  measured  responses. 


The  first  step  in  perfonning  System  Identification  and  model  calculation  is  to  inject  the 
independent  white  noise  signals  into  all  of  the  control  actuators.  This  will  generate  an  /-sample 
input/output  time  history  of  measured  responses,  y,  and  generated  commands,  u.  These  samples 
are  arrayed  into  the  input/output  data  matrix,  V. 
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Given  this  matrix,  the  observer  Markov  parameters,  Y,  can  be  calculated  with  the  following 
equation: 


y  =  YV 


Yt- o' 
.Yt-i. 
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1 

£ 

yt-i+p  +  1- 
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Solve  for  Y  by  applying  the  psuedoinverse  to  V.  The  columns  of  Fean  be  rearranged  to  provide 
the  multi-step  output  predictor. 

Y  =  [t  ft  a±  /?2  a2  ...  Pp  ap\ 

The  resulting  vehicle  model  is  a  weighted  sum  of  the  time  history  data.  This  can  be  used  to 
predict  the  future  responses  given  the  past  history  and  future  commands. 

[yp]  =  x[up]  +  a[yp]  +  P[up] 

yt+o]  ut+o]  Tt-r  rut-i 

:  =  x  :  +  a  :  +  (3  : 

Yt+p.  Mt+p .  yt—p.  Mt-p 

The  next  predicted  response  is  the  first  row  of  the  previous  equation 

Vt+ o  =  Tut+0  +  «iyt-i  +  ftut-i  +  •••  +  «Pyt-P  +  PpU-t-p 

As  with  LQR,  a  control  law  can  be  designed  based  on  an  objective  function, 

n  O0 

J  —  I  (yTQy  +  uTRu)dt 

Jo 

Substituting  the  expanded  equation  for  [y]  yields 

r<  oo  'ut+ o'  Tt-r  ut-i'  ~ut+ o'  yt-r  ~ut-i 

J  =  (  t  :  +a  :  +P  :  ?  t  :  +a  :  +p  :  +  uTRu  )9t 

°  ut+p.  yt—p.  Mt—p.  ut+p.  yt—p.  yt—p. 

Finding  the  local  minimum  of  the  objective  function,  J,  yields  an  optimal  solution  for  the  control 
commands  that  minimize  the  responses.  Taking  the  partial  derivative  with  respect  to  command, 
[u],  and  solving  for  control  commands,  yields  the  following: 

ut+oi  ryt-1!  rut-i 

:  =  —  (xtRt  +  Q )~1xTR  a  '■  +  p  : 

Ut+p.  .Yt-p.  Mt—p 

Defining  the  following  control  law  matrices, 

ac  —  —(ttRt  +  Q)-1x:r/?a, 

Pc  =  —(ttRt  +  Q)-1TrR/?, 

yields  control  laws  that  calculate  the  next  p  commands, 

^t+oi  ryt-n  rut-i 

:  =  «c  :  +  Pc  : 

Ut+p.  yt—p.  Mt-p 

where  ut+p  corresponds  to  the  control  commands  p  steps  into  the  future. 

However,  only  the  next  control  commands  are  of  interest,  thus  the  rows  of  [z<]  corresponding  to 
future  data  can  be  dropped. 
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■yt-i 

~Ut~  1 

[ut+o]  =  ac 

+  Pc 

yt-p. 

Mt-p. 

4.4  GPC  with  Explicit  Disturbance 

Generalized  Predictive  Controller  with  Explicit  Disturbance,  as  described  in  a  subsequent 
publication  by  Kvatemik,  et  al.  (6),  is  similar  to  the  previous  version  of  GPC,  except  that  now 
the  external  disturbances  are  explicitly  modeled.  Time  histories  of  controls  commands, 
measurement  responses,  and,  now,  disturbance  measurements  are  used  to  generate  the  ARX 
model  for  the  system. 

Y  —  [y  t  ^  a1  S2  /?2  a2  ■■■  Sp  Pp  ap\ 

These  are  used  to  predict  future  responses  based  on  past  responses  (including  external 
disturbances)  and  control  commands.  In  this  formulation,  it  is  assumed  that  the  current 
disturbance  can  be  measured. 


[yP]  =  t[up]  +  a[yp]  +  p[up]  +  8[dp] 


yt+ o' 

ut+ o' 

Vt- 1 

ut- 1 

dt-i 

=  T 

+  a 

+  P 

+  6 

yt+p. 

yt+p. 

yt-P. 

Mt—p. 

.d-t-p. 

As  before,  the  next  control  commands  are  calculated  by  minimizing  the  objective  function;  the 
rows  corresponding  to  future  control  commands  are  dropped.  Using  this  extension,  the  previous 
measured  disturbance  infonnation  can  be  used. 


5.  Road  Profiles 


Three  road  profiles  were  used  to  evaluate  the  algorithms:  Sine,  Curb,  and  Perryman  3.  In  each 
case,  the  road  profiles  were  provided  as  an  array  of  ground  distance  and  vertical  heights.  Vehicle 
speed  was  integrated  to  track  the  position  along  the  ground  and  calculate  road  height  at  the  next 
time  step. 

5.1  Sine 

This  road  profile  was  defined  by  a  sine  wave  with  a  half-amplitude  of  20  cm  with  a  wavelength 
of  1  m.  The  intent  was  to  provide  a  single  frequency  disturbance  that  relates  to  vehicle  speed,  in 
this  case,  1  m/s  (3.6  kph)  of  speed  corresponds  to  1  Hz  excitation.  Figure  2  shows  a  portion  of 
this  type  of  road  profile. 
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Sine  Road  Profile 


4  5  6 

Road  Distance  -  m 


Figure  2.  Sine  road  profile. 

5.2  Curb 

This  road  profile  is  defined  by  a  10  cm  step  rise  followed  50  by  meters  of  level  surface,  and  then 
a  10  cm  step  drop.  The  intent  was  to  provide  a  step  input  road  profile  into  the  system.  Note  that 
the  input  was  modeled  as  a  step  input  to  the  tire.  It  does  not  accurately  take  the  geometry  of  a 
physical  tire  rolling  over  curb  into  account,  which  would  effectively  low-pass  filter  the  curb 
input.  Figure  3  shows  this  road  profile. 


Curb  Road  Profile 


Road  Distance  -  m 


Figure  3.  Curb  road  profile. 

5.3  Perryman3 

Perryman  3  is  an  off-road  test  track  located  at  Aberdeen  Proving  Ground.  This  road  profile  is 
intended  to  evaluate  the  algorithms  against  a  challenging,  “real-world”  course.  Figure  4  shows 
this  road  profile. 
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05 


Perryman  3  Road  Profile 


Figure  4.  Perryman  3  road  profile. 


6.  Simulation  Results 


The  vehicle  models  and  controllers  were  implemented  in  Mathworks  Matlab  and  Simulink.  The 
simulation  was  run  using  discrete  time  steps  at  0.01  s  (100  Hz)  using  the  ODE5  (Dormand- 
Prince)  integration  method.  The  simulations  were  run  for  12  s  until  the  model  reached  steady 
state.  The  open  loop  model  was  run  in  parallel  with  closed-loop  LQR,  closed-loop  LQR  with 
Preview  Control,  GPC  with  hnplicit  Disturbances,  and  GPC  with  Explicit  Disturbances  for 
performance  comparison  studies.  The  Preview  Control  block  was  pre-populated  with  the  tPreview 
road  profile  samples  at  the  start  of  the  simulation. 

Initial  conditions  for  the  suspension  were  equal  to  zero  and  the  effects  of  gravity  were  ignored 
such  that  the  system  would  be  at  steady  state  prior  to  the  addition  of  external  disturbances. 
Vehicle  speed  was  held  at  a  constant  value  for  the  entire  run. 

The  following  objective  function  was  used  as  the  metric  for  comparing  algorithm  performance 

]—  low  pass(jyTy  +  uTu  ) 

The  values  for  each  simulation  were  averaged  over  the  12  s  and  plotted  at  multiple  vehicle 
speeds.  The  measured  responses  are  normalized  against  the  open-loop  responses  at  each  speed. 

Normalized  Response  —  1  —  - 'ala° 1  lthm 

J open 

The  actuator  forces  and  powers  are  normalized  against  the  baseline  algorithm,  LQR,  R=l. 

Normalized  Actuator  —  ^ ala°rithm 

Jlqr,r=i 
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This  normalization  allows  a  comparison  of  the  algorithms  despite  the  variation  in  actuator  force 
or  power. 


Normalized  Effectiveness  — 


Normalized  Response  aigorithm 
Normalized  Actuatoralgorithm 


6.1  Algorithms 

A  short  comparison  of  the  four  control  laws  operating  on  two  sets  of  quarter-car  parameters  and 
three  road  courses  (disturbance  inputs)  is  presented  in  the  discussion  below. 


Simulation  data  based  on  the  Light  Vehicle  parameters  is  presented  in  figures  5,  7,  and  12.  These 
figures  were  generated  using  the  following  algorithm  parameters,  as  shown  in  table  3. 


Comparison  of  Algorithms 
Light  Vehicle,  Sine  Course 


Objective  Function 


Normalized  Effectiveness 


Figure  5.  Comparison  of  control  algorithms  for  Light  Vehicle  parameters  on  Sine  course. 
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Table  3.  Algorithm  tuning  parameters  used  for  Light  Vehicle  parameters. 


LQR 

LQR  with  Preview 
Control 

GPC  with  Implicit 
Disturbance 

GPC  with  Explicit 
Disturbance 

sensors 

4  (internal  states) 

4  (internal  states) 

3 

3 

actuators 

Q 

le8*[l  0  1  1] 

le8*[l  0  1  1] 

le5*[l  1  1] 

le5*[l  1  1] 

R 

1 

1 

1 

1 

Preview  Horizon 

- 

20  points 

- 

- 

ARX  model  size 

- 

- 

6 

6 

ID  sample  size 

- 

- 

400  points 

400  points 

ID  disturbance 

- 

- 

10  Hz  sine 

white  noise 

Actuator  Limit 

1.5e4  N 

Simulation  data  using  the  Heavy  Vehicle  parameters  is  presented  in  figures  6,  8,  and  13.  These 
figures  were  generated  using  the  following  algorithm  parameters,  as  shown  in  table  4. 


Comparison  of  Algorithms 
Heavy  Vehicle,  Sine  Course 


Objective  Function 


Figure  6.  Comparison  of  control  algorithms  for  Heavy  Vehicle  parameters  on  Sine  course. 
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Table  4.  Algorithm  tuning  parameters  used  for  Heavy  Vehicle  parameters. 


LQR 

LQR  with  Preview 
Control 

GPC  with  Implicit 
Disturbance 

GPC  with  Explicit 
Disturbance 

sensors 

4  (internal  states) 

4  (internal  states) 

3 

3 

actuators 

Q 

lel2*[l  0  1  1] 

lel2*[l  0  1  1] 

5e9*[l  1  1] 

5e9*[l  1  1] 

R 

1 

1 

1 

1 

Preview  Horizon 

- 

20  points 

- 

- 

ARX  model  size 

- 

- 

6 

6 

ID  sample  size 

- 

- 

400  points 

400  points 

ID  disturbance 

- 

- 

10  Hz  sine 

white  noise 

Actuator  Limit 

1.5e6  N 

Figures  5  and  6  correspond  to  the  Sine  road  profde,  which  represents  a  single  frequency 
disturbance.  In  the  regions  where  the  control  laws  are  valid,  both  GPC  algorithms  show  higher 
perfonnance.  Note  that  the  LQR  based  algorithms  actually  perform  worse  than  the  open  loop 
case  near  vehicle  resonance.  Knowledge  of  disturbance  extends  the  region  over  which  the  control 
laws  are  valid  and  effective,  as  demonstrated  by  GPC  with  Explicit  Disturbances  and  LQR  with 
Preview  Control. 


The  performance  of  the  algorithms  for  the  Curb  road  profile  is  shown  in  figures  7  and  8. 


Comparison  of  Algorithms 
Light  Vehicle,  1 0cm  Curb  Course 


Objective  Function 


- 

- Open 

- LQR 

. LQR  with  20pt  Preview 

- GPC-Explicit  Disturbance 

0  40  50 

Vehicle  Speed  -  kph 
Objective  Function  -  %  Improvement 


Vehicle  Speed  -  kph 
Actuator  Force 


Vehicle  Speed  -  kph 
Normalized  Effectiveness 


Vehicle  Speed  -  kph 


Figure  7.  Comparison  of  control  algorithms  for  Light  Vehicle  parameters  on  Curb  course. 
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Comparison  of  Algorithms 
Heavy  Vehicle,  1 0cm  Curb  Course 

Objective  Function 


10 

8 

6 

4 


- Open 

—  —  LQR 

. LQR  with  20pt  Preview 

- GPC-Implicit  Disturbance 

—  GPC-Explicit  Disturbance 


2  - 

ot 


L 


0  40  50 

Vehicle  Speed  -  kph 
Objective  Function  -  %  Improvement 


_ [ 

80 


Vehicle  Speed  -  kph 


E 

CD 

>  0.5 


Vehicle  Speed  -  kph 


Figure  8.  Comparison  of  control  algorithms  for  Heavy  Vehicle  parameters  on  Curb  course. 

GPC  with  Implicit  Disturbance  performs  much  worse  than  the  open  loop  case  as  shown  in  figure 
9.  This  is  caused  by  the  fact  that  the  measured  value  for  tire  height  has  a  non-zero  offset  relative 
to  the  starting  reference  height. 
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Time  History 
GPC-Implicit  Disturbance 
Light  Vehicle,  10cm  Curb  Course,  18  kph 


Chassis  Accelera  ion  -m/s2 


Chassis  height  -  m 


Time  -  s 


Tire  height  -  m 


Time  -  s 

Actuator  Force  -  N 


o - *r~ 

l/'w 


3 

Time  -  s 


Figure  9.  Time  history  of  GPC  with  implicit  disturbance  for  Light  Vehicle  parameters  on  Curb  course. 

LQR  with  Preview  Control  also  appears  to  have  worse  performance  than  open  loop.  However,  it 
is  clear  from  figure  10  that  it  has  a  much  lower  settling  time.  Also,  in  the  case  of  LQR  with 
Preview  Control,  special  note  must  be  given  to  the  fact  that  the  actuator  commands  begin  well 
ahead  of  the  disturbance,  thus  reducing  accelerations  on  the  vehicle  when  encountering  the  step 
input,  as  shown  in  figure  11.  It  is  inferred  from  results  that  vehicle  speed  is  not  a  major  factor  in 
the  algorithm  performance  since  the  road  input  ignores  the  effect  of  tire  geometry  rolling  up  and 
over  the  curb. 
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Time  History 
LQR  wi  h  20pt  Preview 
Light  Vehicle,  10cm  Curb  Course,  18  kph 


Chassis  Accelera  ion  -m/s2 


Chassis  height  -  m 


Time  -  s 


Tire  height  -  m 


Actuator  Force  -  N 


3 

Time  -  s 


Figure  10.  Time  history  of  LQR  with  preview  control  for  Light  Vehicle  parameters  on  Curb  course. 
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Time  History  (Zoomed) 

LQR  wi  h  20pt  Preview 
Light  Vehicle,  1 0cm  Curb  Course,  1 8  kph 


Chassis  Accelera  ion  -m/s2 


Chassis  height  -  m 

0.5  r  i 

i 


-0  5  L.w _ r _ _ „ _ L_ _ 1 _ 1 _ ft  „  '$ _ I _ Jk _ t _ ,i 

0  0.1  0.2  0  3  0.4  0.5  0.6  0.7  0.8  0.9  1 

Time  -  s 

Tire  height  -  m 


Actuator  Force  -  N 


Time  -  s 

Figure  11.  Time  history  of  LQR  with  preview  control,  close  up  of  step  input. 

The  results  from  the  final  road  profile,  Perryman  3,  are  shown  in  figures  12  and  13.  Again,  GPC 
with  hnplicit  Disturbance  responds  poorly  to  the  input.  LQR,  LQR  with  Preview  Control,  and 
GPC  with  Explicit  Disturbance  perform  well  on  this  aperiodic  road  profile.  Although  the 
Normalized  Effectiveness  appears  similar  for  these  three  algorithms,  it  is  interesting  to  observe 
that  the  objective  function  for  GPC  with  Explicit  Disturbance  is  relatively  flat  over  the  entire 
frequency  range.  This  suggests  that  the  closed-loop  vehicle  response  will  be  independent  of 
vehicle  speed. 
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Comparison  of  Algorithms 
Light  Vehicle,  Perryman  3  Course 


Objective  Function 


Figure  12.  Comparison  of  control  algorithms  for  Light  Vehicle  parameters  on  Perryman  3  course. 
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Comparison  of  Algorithms 
Heavy  Vehicle,  Perryman  3  Course 


Objective  Function 


Figure  13.  Comparison  of  control  algorithms  for  Heavy  Vehicle  parameters  on  Perryman  3  course. 

The  quarter-car  model  assumed  that  the  internal  states  were  available  for  feedback  control. 
Several  of  these  states  are  not  good  choices  for  optimal  control  or  practical  implementation.  For 
instance,  tire  height  above  the  starting  reference  point  skews  the  objective  function  for  a  step 
input  and  chassis  vertical  velocity  is  difficult  to  measure. 

6.2  Preview  Horizon 

Preview  control  requires  knowledge  of  future  disturbances.  These  future  disturbances  can  be 
measured  or  modeled  (6).  Figure  14  shows  the  effects  of  increasing  the  preview  horizon  from  0 
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(basic  LQR  algorithm)  to  100  points.  Very  few  preview  points  are  required  to  achieve  large 
reductions  in  the  objective  function;  large  numbers  of  preview  points  do  not  improve  the 
perfonnance.  This  corresponds  well  with  the  results  reported  by  El  Madany  (2). 


Comparison  of  Algorithms 
Heavy  Vehicle,  Perryman  3  Course 


Objective  Function 


Figure  14.  Comparison  of  effect  of  preview  horizon  on  LQR. 

Practical  sensors  to  measure  this  value  do  not  exist  outside  of  a  controlled  laboratory.  No  effort 
was  made  to  identify  and  model  a  sensor  to  directly  gather  disturbance  data  or  preview  data.  All 
measurements  assumed  no  noise  or  offset -bias. 
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6.3  System  Identification 

The  LQR  formulation  requires  complete  and  perfect  knowledge  of  the  vehicle  model  and  internal 
states  prior  to  calculation  of  the  control  laws.  While  this  is  not  practical  for  most  applications,  it 
does  provide  a  common  algorithm  to  compare  during  simulation.  This  limitation  can  be 
mitigated  with  the  Linear  Quadratic  Gaussian  (LQG)  fonnulation,  which  includes  a  Kalman 
filter  to  adjust  an  imperfect  vehicle  model.  This  is  not  considered  in  the  scope  of  this  currently 
reported  research  effort. 

Generalized  Predictive  Control  algorithm  relies  on  a  system  identification  to  generate  control 
laws.  Typically,  a  given  system  identification  is  valid  over  a  small  frequency  range  about  the 
primary  excitation  frequency  of  the  disturbance.  The  robustness  of  the  system  identification  was 
evaluated  for  different  types  of  model  disturbances.  Figures  15  and  16  show  the  effects  of 
system  identifications  under  different  disturbance  types. 

As  shown  in  figure  15,  GPC  with  Implicit  Disturbances  generates  control  laws  that  are  valid  over 
a  range  of  20  Hz  and  the  valid  range  shifts  with  the  System  Identification  disturbance  frequency. 
Surprisingly,  when  GPC  with  Implicit  Disturbances  generates  control  laws  with  disturbances 
composed  of  multiple  frequencies  (either  white  “noise”  or  a  signal  composed  of  1,  5,  7,  11,  and 
17  Hertz),  the  control  laws  are  not  valid  over  the  entire  frequency  range.  This  suggests  that  new 
system  identifications  are  required  at  higher  vehicle  speeds. 
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Evaluation  of  Disturbance  frequency  content  on  System  Iden  ification 
GPC  with  Implicit  Disturbance 
Light  Vehicle,  course:  Sine  Course 


Objective  Function 


100 


Objective  Function  -  %  Improvement 


Normalized  Effectiveness 


Figure  15.  Comparison  of  effect  of  frequency  content  of  disturbance  for  GPC  with  implicit  disturbances 
system  identification. 


The  performance  of  GPC  with  Explicit  Disturbance,  as  shown  in  figure  16,  improves  with  the 
addition  of  white  “noise”  or  composite  disturbances.  In  each  case,  it  is  assumed  that  the  actuator 
excitation  and  disturbances  are  completely  known  (no  noise). 
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Evaluation  of  Disturbance  frequency  content  on  System  Iden  ification 
GPC  with  Explicit  Disturbance 
Light  Vehicle,  course:  Sine  Course 


Objective  Function 


Actuator  Force 


Normalized  Effectiveness 


Figure  16.  Comparison  of  effect  of  frequency  content  of  disturbance  for  GPC  with  explicit  disturbances 
system  identification. 

The  weakness  of  GPC  with  hnplicit  Disturbances  can  be  overcome  by  performing  real-time 
system  identification  at  new  operating  speed.  A  new  system  identification  can  be  performed 
while  the  control  system  is  already  in  closed-loop.  This  allows  the  algorithm  to  use  a  control  law 
optimized  over  the  current  operating  range.  This  approach  will  outperform  GPC  with  Explicit 
Disturbances  with  a  given  periodic  disturbance  as  shown  in  figure  17. 
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Multiple  System  Identification  at  Speed 
Light  Vehicle,  Sine  Course 


Objective  Function 


Objective  Function  -  %  Improvement 


Actuator  Force 


Normalized  Effec  iveness 


Figure  17.  Effects  of  multiple  system  identifications  for  GPC  with  implicit  disturbance  of  closed  loop 
system  in  real-time. 

6.4  Control  Weighting 

As  shown  in  figures  18,  19,  20,  and  21,  when  the  improvement  of  the  objective  function  is 
normalized  by  the  actuator  force,  the  variation  of  the  weighting  factors,  (Q  and  R),  does  not 
greatly  affect  the  performance  of  an  algorithm.  This  property  is  used  to  compare  the 
effectiveness  of  the  different  algorithms  despite  variations  in  actuator  forces,  sensor  magnitudes, 
and  weighting  factors.  This  statement  is  not  valid  when  the  control  laws  are  unable  to  control  the 
system. 
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Improvement/Norm  Force 


Comparison  of  Control  Weights 
LQR 


Light  Vehicle,  Sine  Course 


Objective  Function 


Objective  Function  -  %  Improvement 


Actuator  Force 


Normalized  Effectiveness 


Figure  18.  LQR  algorithm  at  various  control  weights. 
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Improvement/Norm  Force 


Comparison  of  Control  Weights 
LQR  wi  h  20pt  Preview 
Light  Vehicle,  Sine  Course 


Objective  Func  ion 


Objective  Func  ion  -  %  Improvement 
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Figure  19.  LQR  with  20  point  preview  control  at  various  control  weights. 
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Improvement/Norm  Force 


Comparison  of  Control  Weights 
GPC-Implicit  Disturbance 
Light  Vehicle,  Sine  Course 


Objective  Func  ion 


Objective  Function  -  %  Improvement 


Actuator  Force 


Normalized  Effectiveness 


Figure  20.  GPC  with  implicit  disturbances  at  various  control  weights. 
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Comparison  of  Control  Weights 
GPC-Explicit  Disturbance 
Light  Vehicle,  Sine  Course 
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Figure  21 .  GPC  with  explicit  disturbances  at  various  control  weights. 


7.  Conclusions 


Knowledge  of  the  disturbance  clearly  shows  an  improvement  in  the  performance  of  the  control 
laws.  Both  LQR  with  Preview  Control  and  GPC  with  Explicit  Disturbances  outperfonn  the 
versions  of  the  algorithms  without  disturbance  knowledge.  A  relatively  short  distance  of  preview 
disturbance  is  required  to  show  improvements.  Longer  distance  preview  distances  do  not 
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improve  performance  much.  It  is  shown  that,  for  a  given  algorithm,  the  control  weights  reduce 
sensor  objective  function  at  the  same  rate  that  it  increases  actuator  objective  function.  Therefore, 
control  weights  can  be  normalized  for  comparison  across  algorithms.  Though  control  laws  based 
on  GPC  with  Implicit  Disturbance  were  valid  over  a  very  narrow  speed  range;  the  system  can  be 
re-identified  periodically  to  achieve  the  best  performance  of  all  of  the  algorithms. 


8.  Future  Work 


Generalized  Predictive  Control  with  Explicit  Disturbances  will  be  extended  to  include  Preview 
Control,  and  non-linear  system  identification  will  be  incorporated  to  allow  for  semi-active 
suspension  system  development.  These  controllers  will  be  implemented  on  high  fidelity  multi¬ 
body  dynamics  models  and  will  be  experimentally  evaluated  on  a  physical  suspension  test  rig. 

A  two-degree-of-freedom  suspension  test  rig  has  been  fabricated,  as  shown  in  figure  22.  This 
will  be  used  to  validate  aspects  of  the  control  laws  in  hardware-in-the-loop  control.  This  test  rig 
has  variable  sprung  and  unsprung  masses  of  up  to  250  kg,  representing  a  1000  kg  vehicle.  The 
tire  and  suspension  components  can  be  reconfigured  using  commercial-off-the-shelf  components 
or  custom  components  as  needed. 


Actuator  Loadcell 


%  Total  Chassis  Mass 
(15-250  kg)' 


COTS  or  Custom 
Suspension  Components 


LVDT 


Tire  Mass 
(5-150  kg)' 


Tire  Model 


Variable  Tire  Stiffness 


LVDT 


Road  Actuator 
±75  cm  (a)  2  Hz- 
±7  cm  @  20  Hz 


Tire  Loadcell 


Figure  22.  2  DOF  suspension  test  rig. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


2  DOF 
A 

accent  - 

accent  - 

accent • 

accent  •• 

ARE 

ARL 

ARX 

B 

C 

d 

F 

GPC 

J 

K 

K 

LQG 

LQR 

m 

MR 

N 

P 

Q 

R 

subscript  c 
subscript  d 


two-degree-of-freedom 

state-space  state  matrix 

approximated  value 

first  row  of  Markov  coefficient  matrix 

first  derivative  with  respect  to  time 

first  derivative  with  respect  to  time 

Algebraic  Riccati  Equation 

U.S.  Anny  Research  Laboratory 

Autoregressive  with  Exogenenous 

state-space  control  input  matrix 

damper  coefficient 

external  disturbance  vector 

external  control  force 

Generalized  Predictive  Control 

objective  function 

spring  Constant 

Feedback  Gain  Matrix 

Linear  Quadratic  Gaussian 

Linear  Quadratic  Regulator 

mass 

magneto-rheological 
cross-coupled  weighting  matrix 
solution  to  Algebraic  Riccati  Equation 
response  weighting  matrix 
control  command  weighting  matrix 
property  of  control  command 
property  of  external  disturbance 
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subscript  S 
subscript  t 
subscript  U 
superscript  T 
t 
u 

VTD 

x 

y 

z 

a 

P 

5 

x 


property  of  sprung  or  chassis  mass 
property  of  current  time  step 
property  of  unsprung  or  tire  mass 
matrix  transpose 
current  time  step 
control  command  vector 
Vehicle  Technology  Directorate 
state  vector 

measured  response  vector 
element  of  internal  state  vector 
Markov  response  coefficient  matrix 
Markov  control  coefficient  matrix 
Markov  disturbance  coefficient  matrix 
Markov  future  command  coefficient  matrix 
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NO.  OF 

COPIES  ORGANIZATION 

1  ADMNSTR 

ELEC  DEFNS  TECHL  INFO  CTR 
ATTN  DTICOCP 

8725  JOHN  J  KINGMAN  RD  STE  0944 
FT  BEL  VOIR  VA  22060-6218 

1  CD  OFC  OF  THE  SECY  OF  DEFNS 
ATTN  ODDRE  (R&AT) 

THE  PENTAGON 
WASHINGTON  DC  20301-3080 

1  US  ARMY  RSRCH  DEV  AND  ENGRG 

CMND 

ARMAMENT  RSRCH  DEV  &  ENGRG 
CTR 

ARMAMENT  ENGRG  &  TECHNLGY 
CTR 

ATTN  AMSRD  AAR  AEF  T 
J  MATTS 
BLDG  305 

ABERDEEN  PROVING  GROUND  MD 
21005-5001 

1  PM  TIMS,  PROFILER  (MMS-P) 
AN/TMQ-52 
ATTN  B  GRIFFIES 
BUILDING  563 
FT  MONMOUTH  NJ  07703 

1  US  ARMY  INFO  SYS  ENGRG  CMND 

ATTN  AMSEL  IE  TD  A  RIVERA 
FT  HUACHUCA  AZ  85613-5300 

1  COMMANDER 

US  ARMY  RDECOM 
ATTN  AMSRD  AMR 
WC  MCCORKLE 
5400  FOWLER  RD 

REDSTONE  ARSENAL  AL  35898-5000 

1  US  GOVERNMENT  PRINT  OFF 

DEPOSITORY  RECEIVING  SECTION 
ATTN  MAIL  STOP  ID  AD  J  TATE 
732  NORTH  CAPITOL  ST  NW 
WASHINGTON  DC  20402 


NO.  OF 

COPIES  ORGANIZATION 

7  US  ARMY  RSRCH  LAB 

ATTN  RDRLVTA  HL  EDGE 
ATTN  RDRLVTA  J  PUSEY 
ATTN  RDRLVTM  D  LE 
ATTN  RDRLVTM  M  MURUGAN 
ATTN  RDRL  VTM  R  BROWN 
(3  HCS) 

ABERDEEN  PROVING  GROUND  MD 
21005 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRL  CIM  G  T  LANDFRIED 
BLDG  4600 

ABERDEEN  PROVING  GROUND  MD 
21005-5066 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRLVTA  JA  BORNSTEIN 
ABERDEEN  PROVING  GROUND  MD 
21005-5066 

3  US  ARMY  RSRCH  LAB 
ATTN  IMNEALCHRR 
MAIL  &  RECORDS  MGMT 
ATTN  RDRL  CIM  L  TECHL  LIB 
ATTN  RDRL  CIM  P  TECHL  PUB 
ADELPHI  MD  20783-1197 

TOTAL:  19(1  ELEC,  1  CD,  17  HCS) 


33 


Intentionally  Left  Blank. 


34 


